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Bull. Inst. Zool., Academia Sinica 28(3): 


183-189. The 


mating behavior in flies for a closed as well as an open environment are modeled 


stochastically through a Markov chain representation. 


The stochastic approach 


provides not only the mean or macroscopic information about the phenomenon under 


consideration but also an estimation of its fluctuating characteristic. 


The applica- 


bility of the present stochastic model is justified by analyzing the available experi- 


mental data. 
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A and increasing interest in 
modeling the mating behavior in flies 
has been observed (Eckstrand and Seiger, 
1975; Taylor, 1975; Kence and Bryant, 
1978; Dowse et al., 1984; Dowse et al., 1986; 
Hsu and Tseng, 1988). Most of the 
mating behaviors, however, are described 
in a deterministic manner (see, Taylor, 
1975; Dowse et al., 1986; Hsu and Tseng, 
1988). That is, the behavior of a dynamic 
system is represented by its mean or 
macroscopic characteristic. Thus, given 
the initial condition of the system, its 
status at an arbitrary point on the time 
scale is uniquely determined. In contrast 
to a deterministic model, a dynamic 
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system is portrayed through probabilistic 
statements in a stochastic model. The 
solution to such a model yields probability 
distributions describing the variations of 
the random variables representing the 
system. In general, the mean or the 
first moment of the probability distri- 
bution of the random variable of a 
stochastic model reduces to the result 
predicted by the corresponding deter- 
ministic model. Therefore the former is 
capable of providing more detailed in- 
formations about a dynamic system than 
the latter. In other words, a stochastic 
model is a generalization of the cor- 
responding deterministic model. 

Unlike a lifeless entity, each individual 
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fly possesses its characteristic. A deter- 
ministic representation is unable to 
reflect this fact since flies are viewed 
as identical entities having the mean 
characteristic of the population of the 
system. On the other hand, the pro- 
babilistic approach in a stochastic re- 
presentation provides a possibility of 
portraying the random nature inherent 
in the phenomenon. Furthermore, it is 
well known that the mean and the 
variance of a random variable is on the 
same order of magnitude. Therefore if 
the number of flies in a dynamic system 
is not large enough, the mean of the 
random variable representing the number 
of flies is not suflicient in describing its 
behavior. The random nature of the 
phenomenon is especially significant in 
this case. In the present study, we 
extend our previous deterministic analysis 
- (Hsu and Tseng, 1988) to the correspond- 
ing stochastic counterpart, in the form 


of a Markov chain (Bharucha-Reid, 1960). - 


An attempt is made such that the degree 
of the fluctuation in the number of mated 
fly pairs can be estimated. 


MARKOV CHAIN REPRE- 
SENTATION 


The transition stage mating me- 
chanism proposed by Hsu and Tseng 


1 
1| pu(m, m+l) 
Pim, m+1)= 2 | palm, m+1) 
3 L bam, m-+1) 
The rate of formation of the fly pairs in 
the transition stage is assumed to be 
proportional to the product of the 


number of males and the number of 
females; and the rate of formation of 
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(1988) is adopted to describe the mating 
behavior in flies. According to this 
mechanism, a male fly and a female fly, 
when encountered, need to pass through 
a transition stage before they become a 
mated pair. Symbolically, we have 


M+4-F-=> MF*— MF (1) 


where M, F, MF*, and MF denote, re- 
spectively, male flies, female flies, fly 
pairs in the transition stage, and mated 
fly pairs. It is defined that flies in states 
1, 2, and 3, represent male flies, fly pairs 
in the transition stage, and mated fly 
pairs respectively. The number of males, 
pairs in the transition stage, and mated 
pairs, are denoted by nı, na and 13, 
respectively. The number of females, nr, ' 
can be evaluated by the following ex- 

pression: bp 


np=F,—(M,—n1) 
=F—M+n1 ( 2) 


where M, and F, are the males and 
females initially present in the system. 
respectively. Define p;,(m, m+1) as the pro- 
bability that a fly in state i at time mAt 
will be in state j at time (m+1)At where 
At is a unit time interval. The evolution 
of the mating phenomenon is thus 
described by the following transition 
probability matrix P(m, m+1): 


2 3 
Die(m, m+1) pulm, m+1) 
belm, m+1) palm, m+1) (3) 
Da(m, m+1) Dss(m, m+1) 


the number of mated pairs is proportional 
to the number of fly pairs in the transi- 
tion stage (Hsu and Tseng, 1988). These 
assumptions lead to the following transi- 
tion probability matrix: 


2 3 
kınr(m)At 0 | 
1—k, At hk, At (4) 
0 - Jj 


MATING BEHAVIOR IN FLIES 


where k, and & are constant. Clearly, 
if there are n,(m) flies in state 7 at time 
mât, these flies will be in one of the 
states j at time (m-+1)Af with probability 
Din, m+l), i, 71, 2, 3. That is, ni(m+1) 
is multinomially distributed among the 
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states 1, 2, and 3. Consequently, the 
conditional mean and the variance of a; 
at (m+1)At, given the value of n; at time 
mAt, i, ¡=1,2,3, are respectively (Rohatgi, 
1976), 


3 
Eln;(m-+1) |nı(m), i=1, 2, 3] = 2i n(m) pulm, m+1), j=1, 2, 3. (5) 
and | 
Var[nj(m+1)|u:(m), i=1, 2, 3] 
=Y n;(m) Di On, m+1)[1—p;;(m, m+1)1, j=l, l, 2, 3. (6) 
the following differential equations: 
DATA ANALYSIS 
UE a, (10). 
The applicability of the present dt 
stochastic model is examined through dn, _, Lk (11) 
analyzing the available experimental ae 
results. For illustration, the data reported dns | 
by Dowse et al. (1986) are adopted. The = Ren (12) 


mating experiment was conducted with 
one strain of female and one strain of 
male Drosophila flies. Their experimental 
data along with the values predicted by 
the present stochastic model are pre- 
sented in Figs. 1 through 3. Also shown 
in these figures is an approximate 95% 
confidence interval, calculated by E[.]+ 
1.96Var[. 11. The estimation of the 
values of the adjustable parameters kı 
and k, are based on equation (5). The 
least sum of error squares criterion is 
adopted. | i 


DISCUSSION 


Note that the time scale in a Markov 
chain model is discrete. If the mating 
behavior under consideration is expressed 
in a deterministic manner, we have 


nı(m+1) =nı(m) —kını(m) nr(m)At 


(7) 

No(m+1)=n2(m) +h, nı(m) nr(m)At 
— k, nolm)At ( 8 ) 
ns (m+1l)=ns(m)+kın,(m)At (9) 


where np=Fọ— Mtn, In the limit as 
At-0 equations (7) through (9) reduce to 


Thus the deterministic counterpart of 
the present stochastic model is equivalent 
to that proposed by Hsu and Tseng 
(1988). This is also justiied by the 
variation of the mean of the number of 
mated pairs illustrated in Figs. 1 through 
3. The present stochastic model is 
capable of providing more detailed in- 
formation about the dynamic system. 
For instance, ‘the variation of the 
number of flies, which is necessary in 
constructing the confidence bands shown 
in Figs. 1 through 3. It should be 
pointed out that np is a random variable, 
and therefore, its value at an arbitrary 
time. mÁt can not be determined except 
at m=0. Hence in evaluating the values 
of the elements in the transition pro- 
bability matrix defined by equation (4), 
nr is replaced by its mean or expected 
value Elnr]. This procedure has the 
effect of reducing the magnitude of the 
variances of the number of flies. The con- 
fidence interval shown in Figs. 1 through 
3 should thus be viewed as a lower bound, 
i.e. the actual bands can be much wider 
than those presented. 
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Fig. 1. Experimental data of mating of D. simulans with Mo=8 and F.=16 (Dowse et.al., 1986) 
along with the predicted results by the present model, k1==3.688 x 10-3/min, 42=6.6237 X 10-!/min. 
——: mean number of mated pairs; ---: approximated 95% confidence interval. 


The present stochastic model can be 
generalized without much difficuly if the 
‘mating of flies occurs in an open en- 


vironment. Let us consider the case if 
flies are allowed to enter or leave a 
1 
2 0 
l= 
P(m, m+1) 3 0 
di 0 


where ks, k, and k; are constant. A fly 
or fly pair in state 4 denotes it leaves 
the volume. Let Ni(m) be the rate of 
flies in state ¿ entering the volume at 


specified volume in space. Suppose that 
the rate of fly leaving this volume is 
proportional to its number in the volume. 
Then the transition probability matrix 
takes the following form: 


-A 3 4 
kı nr (m)At 0 k At | 
1—[k.+k,]At ko At k, At (13) 
0 1—k At ks At 
0 0 1 
time mAt. The conditional mean and 


the variance of n; at (m+1)At, given the 
values of n; and N; at time mAt, i,j= 
1, 2,3, are, respectively, 
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Fig. 2. Experimental data of mating of wild-type D. melanogaster (a) Mo=8, Fo=32; (b) M= 

32, Fo=8 (Dowse et al., 1986) along with the predicted results by the present model: ki= 

5.2429 x 10-2/min, k2=1.8793X10-*/min in case (a); ki=5.8271X10-?/min, k2=1.0095x10-*/ 

min in case (b). ——: mean number of mated pairs; ---: approximated 95% confidence 


interval. 
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Time (min) 
(a) and (b) show the experimental data of mating of D. simulans and D. melanogaster vg 
males when Mo=Fo=8 (Dowse et al., 1986) along with the predicted results by the present 
model: %1=6.4407x10-*/min, 2=8.6517X10-1/min in case (a); ki=6.2025X10-*/min, k:= 
9.5079 10-!/min in case (b). ——: mean number of mated pairs; ---: approximated 
95% confidence interval. 
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EIln;(m-+1) n: (m), N;(m), i=1, 2, 3] 
= Ynılm) pulm, DEN, j=1, 2 3. (14) 
and 
Var[n;(m-+1)|n;(m), Ni(m), i=1, 2, 3] 
= Si n(m) buon, m+DLL-Buiom m+], j=1, 28 (15) 
where N;(m), i=1,2,3, is assumed to be correlated, i,j=1,2,3, equations (14) and 


deterministic. If N;(m) is also randomly (15) should be modified as 
distributed with n,(m) and Ni(m) un- 


Eln;(m+1) Inı(m), Ni(m), i=1, 2, 3] 


= E num) piso, m+1)+ELN;(m)], j=1, 2 3. (16) 
and 
Var[n;(m+1) |n:(m), Ni(m), ¿=1, 2,3] : 
= Y nlm) buon, m+DLI-bi0m, mW 1+VarlNi(m)] 35-123. A7) 
where EIN;(m)] and Var[N;(m)] denote, density and mating rates in Drosophila pseudo- 
respectively, the mean of N;(m) and the obscura. Evolution 29: 287-295. 


Hsu, J.P. and S.J. Tseng (1988) Kinetic analysis 


variance of N;(m). ur ce 
of mating behavior in flies: a transition stage 


model. J. Theor. Biol. 130: 175-182. 
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